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Abstract 

We consider a 3+1 dimensional field theory at a Lifshitz point for a dy- 
namical critical exponent z=3, with a scalar and a fermion field coupled via a 
Yukawa interaction. Using the non-perturbative Schwinger-Dyson approach 
we calculate quantum corrections to the effective action. We demonstrate 
that a first order derivative kinetic term as well as a mass term for the 
fermion arise dynamically. This signals the restoration of Lorentz symme- 
try in the IR regime of the single fermion model, although for theories with 
more than one fermionic species such a conclusion will require fine-tuning of 
couplings. The limitations of the model and our approach are discussed. 
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1 Introduction 



Quantum field theory models, in which the UV behavior is governed by a 
Lifshitz-type fixed point have attracted attention recently, as their renormal- 
ization properties appear significantly improved, compared to models with a 
Lorentz symmetric Gaussian fixed point. A novel quantum gravity model, 
which claims power counting renormalizability, has been formulated recently 
by Horava in [H |2] . This scenario is based on an anisotropy between space 
and time coordinates, which is expressed via the scalings t — )■ b^t and x — >■ bx, 
where 2 is a dynamical critical exponent. For z 1 the UV behavior of the 
model is governed by a nonstandard Lifshitz fixed point, while for z = 1 we 
recover the well known Gaussian fixed point. Note that in the Horava model, 
z = 3 is chosen. 

Horava gravity has stimulated an extended research on cosmology and 
black hole solutions, see for example [Sj IH [5l El [TJ [8]. We would like to 
note that Horava gravity is a non-relativistic theory, however it is expected 
that general relativity is recovered in the IR limit. Moreover, some possible 
inconsistencies on Horava gravity have been remarked in [9l HOl [H] , but they 
will not be discussed here. 

Independently of general relativity, quantum field theory models in fiat 
space-time with anisotropy have been studied as well. For example, a thor- 
ough study on renormalization properties of models with a Lifshitz-type fixed 
point, is presented in [111 [131 [HI [15] , ci-nd the Standard Model in this Lorentz 
violating approach is examined in [16]. Also, the renormalizability of scalar 
field theory at the Lifshitz point is examined in [17], and in [18] renormaliz- 
able models with a Lifshitz fixed point are constructed, whereas a renormal- 
izable asymptotically free Yang Mills theory, in 4+1 dimensions, is given in 
[I9] . As far as dynamical mass generation is concerned, a four-fermion inter- 
action has been studied in the framework of Lifshitz-like theories pO] , where 
the authors find a gap equation for the fermion mass, and the CP^~^ model 
at the Lifshitz point is discussed in [21]. In addition, [22] shows some per- 
turbative properties of Lifshitz-like theories containing scalars and fermions, 
where an extension of supersymmetry to a Lorentz non-invariant theory is 
studied. For a presentation of renormalization group equations in the case of 
a scalar field, see [23]. Finally, in [21] a U(l) Gauge theory in 2+1 dimensions 
with 2; = 2 is considered. 

The literature mainly deals with perturbative studies, and our aim here is 
to study a simple field theory model, in the framework of the non-perturbative 
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Schwinger-Dyson approach. In particular, we consider a Lifshitz-type model, 
in flat space time and in 3+1 dimensions, for a dynamical critical exponent 
z=3, with a scalar and a fermion field interacting via a Yukawa coupling. For 
the construction of the bare action of the model, we use only the quadratic 
marginal operators (kinetic terms), with dimension six, plus a Yukawa in- 
teraction term with a dimensionful coupling. Note that the construction of 
more complicated models, including other marginal and relevant operators 
(for 2; = 3) is possible. However, in this work we will restrict our study to 
a Yukawa interaction only, in order to deal with a tractable system of equa- 
tions, describing the dynamical generation of mass and Lorentz symmetry 
for fermions. 

We stress that, for the specific model we examine, the generation of a term 
of the form Xiipj'^df^tp is enough for the restoration of Lorentz symmetry in 
the IR limit of the fermionic sector. However, as Ref.[23] points out, the issue 
of Lorentz symmetry restoration becomes problematic when one considers 
models with more than one species of particles, and we will discuss this 
point before the conclusion. 

To summarize, this study is based on the Schwinger-Dyson approach, for 
which we derive in Appendix A the corresponding equation for the fermion 
self energy. The latter is parametrized by two dressed parameters, m/ and A, 
via the operators iv?jipip and Xiip'-j^d^ip, and the corresponding self consistent 
equations are solved. Note that the parameter A controls the restoration of 
Lorentz symmetry in the fermionic IR sector. The coupled evolutions of 
these two parameters with the Yukawa coupling is presented in figJJl where 
it is found that there exists a critical value for the coupling, above which 
quantum corrections can generate simultaneously a Lorentz invariant kinetic 
term and a mass for fermions. We also comment on the physical relevance of 
our model and the limits of our approximations in Appendix B and Appendix 
C. 

2 Free systems 

We construct in this section the free scalar and fermion models, and derive 
the corresponding propagators which will be used for the loop calculations 
in the next section. 
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2.1 Scalar field 

Here we remind the reader the construction of an anisotropic scalar model, 
in D + 1 dimensions, starting with the action 

Sb = lj dtd^x - (-A)^ 0) , (1) 

where a dot over a letter represents a time derivative. The action ([T]) describes 
a free scalar theory, with the following mass dimensions 

[^1 = -1 [t] = -^ ['^] = ^' 
and leads to the following equation of motion 

0+(-A)^0 = O. (3) 
We look for a solution by assuming the separation of variable 

0(t,x) =e(t)exp{z p-x}, (4) 

which leads to 

e+(p're = o. (5) 

We obtain 

^ = ^0 exp {±ituj) , u = (p2) 5 , (6) 

where ^0 is a constant, such that the solutions of the form of eq.(jl]) represent 
plane waves in D+1 dimensions, and the Feynman propagator for the scalar 
field, which will be used in order to calculate loop diagrams, is 

Gb{u;, p) = — — . ^, , (7) 

where [u] = z. If we include a mass term — |m^^02 in the action of eq.([T]) 
the scalar field propagator is modified as 



^ — (p ) — + i£ 



where [nih] 
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2.2 Fermionic field 

The action for the free fermionic model is 

Sf = j dtd^x l^nV + ^ (-A)^ (h'^dk) , (9) 

where we have included only quadratic marginal operators which correspond 
to a Lifshitz fixed point at the ultraviolet. A dimensional analysis gives 

[x'] = -i [t] = -z m = ^, (10) 

and the equation of motion is: 

,/^ + (_A)^(z7'^9fc)7/' = 0. (11) 

We make the following ansatz for the solution of the above equation 

ilj{t, x) = 9{t)'il:p exp{ip ■ x}. (12) 

where the spinor part -^p is normalized according to the equation '4'p^p = 1- 
If we multiply with the Hermitian conjugate we obtain 

^ + {-AYij = (13) 

The solution f[T^ should satisfy eq. ffT^ . hence we obtain 

e{t) = Ooexp {±itu) , u = (p2)i (14) 

where Oq is a constant, such that the solutions (fT2|) represent plane waves in 
D+1 dimensions. The Feynman propagator for the fermion field is 

Gf{u,p) = ^ (15) 

UJ^ — (p2) 2 (p . 7) _)_ ^5 

. ^7°-(p')^(p-7) 

— 1 

— (p^)^ + ie 

where [u] = z. We can include the mass term —m^jipil) in the action ([9]), 
where [mj] = 1, as well as an additional quadratic ternjfl Xif) (i'^'^dk) ip, where 
[A] = z — 1, such that the fermion propagator is finally 



(p2) 2 + A (p ■ 7) + 
Gf{u;,p) = i -2 (16) 



z-l 



p2)— + A 



p2 — my + te 



^This term is quadratic in the fermion field, but it is not marginal for z 7^ 1. 
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3 Dynamics 



3.1 Model and Schwinger Dyson equations 

We now consider the simplest interaction between scalars and fermions in the 
Lifshitz context, through a Yukawa couphng, and start with the following 
bare action 

S = J dtd^xS^ijij^ij + ij{-A)'^ {i-f''dk)i^ (17) 

where the coupling constant has dimension [g] = In the framework of 

the gradient expansion, we will consider quantum corrections up to the first 
order in momentum only, such that we will look at the corrections to the 
scalar mass, and will allow the dynamical generation of a fermion mass term 
—m^j-tpip and of the additional first order fermionic kinetic term Xipii'-y^dkjip, 
in order to study the restoration of Lorentz invariance for fermions. 

In the action f|T71) . we start with a bare scalar mass in order to absorb 
the only UV divergence which will appear, as we will see, in the corrections 
to the scalar mass. No UV divergence will appear in the fermion self energy, 
due to the higher order derivatives, and for this reason mj and A can be 
taken equal to zero in the bare action. We note here that, also because of 
higher derivatives, the UV divergence we will find in the corrections to the 
scalar mass is logarithmic for D = z = ?>, and not quadratic as it is in a 
Lorentz-invariant theory. 

We will use here the Schwinger Dyson approach to calculate the fermion 
and scalar self energies, which is non-perturbative and represents a resumma- 
tion of graphs, avoiding IR divergences, because of the presence of a fermion 
mass and first order derivative kinetic term, both generated dynamically. 
Also, studies of dynamical mass generation usually lead to a mass which is 
non-analytical in the coupling constant, which cannot be found with a naive 
loop-expansion, and one therefore needs a non-perturbative approach. We 
show in the Appendix that the corresponding Schwinger-Dyson equation for 
the fermion self energy S/ = — is 

i:f = iggfQG,, (18) 
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Figure 1: The fermion self energy given by the Schwinger Dyson equation in the 
rainbow approximation. A sohd thick hne represents the dressed fermion propaga- 
tor, a sohd thin hne the bare fermion propagator, and a dashed hne represents the 
dressed scalar propagator (which, in our approximation, is like the bare propaga- 
tor, but with the renormalized mass instead of the bare one). As one can see, the 
fermion self energy is obtained as a resummation of an infinite number of graphs, 
which is at the origin of the non-perturbative feature of the results. 

where Qf,Qb and G are respectively the dressed fermion propagator, the 
dressed boson propagator and the dressed vertex. The equation f ITS]) is self 
consistent, since it displays the dressed quantities on both sides, and therefore 
corresponds to a resummation of all quantum corrections (see fig.([T])). 

Using the exact equation ( [T8l) . we can study the dynamical generation 
of mass and first order derivative terms for fermions, and we will make the 
following assumptions: 

• We neglect quantum corrections to the vertex, which corresponds to 
the so-called ladder or rainbow approximation j2Bj, and we therefore 
consider G ~ The corresponding partial resummation provided by 
the Schwinger-Dyson equations f|T8|) is the dominant one for the study 
of dynamical mass generationl^. We show in Appendix B that this ap- 
proximation is well controlled in the regime where we observe dynamical 
mass generation; 

• We also neglect the renormalization of the bare fermion kinetic term, 
which is consistent in the framework of the gradient expansion, if we 

^ There is in principle a infinite tower of Schwinger-Dyson equations, which are self 
consistent equations for every 71-point function, each involving the n + 1-point function. A 
given truncation of this tower of coupled equations consists then is a specific resummation 
of graphs for each correlation function. 
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take into account first order derivative corrections to the fermion dy- 
namics only; 

• Also because of the gradient expansion, we consider a momentum- 
independent dynamical mass, since the latter would be quadratic in 
the momentum. In addition, the dominant contribution of the loop 
integral appearing in the Schwinger- Dyson equation flTSl) arises from 
low momentum, since no UV divergence occurs in the calculation of 
the fermion self energy. This approximation is discussed in Appendix 
C. 

In what follows, we will concentrate on the case D = z = ?>. 



3.2 Scalar sector 

It can be shown, as done in the Appendix for the fermion self energy, that 
the Schwinger Dyson equation for the scalar self energy reads 

= iT{g^^ - G^^} = igiiigfQgj}. (19) 

As we will see in the next subsection, the operators ifjij) and ijj{i'j''dk)ip "will be 
generated dynamically, such that we assume here that the dressed fermion 
propagator has the form ( fT6l) . Qf = Gf, where rrif and A are generated 
dynamically. The scalar mass, after a Wick rotation, is then obtained from 
eq. (fT9|) for vanishing momentum, which reads 

m<^-m^ = Aa^ r - ! ^^ + (P^ + A)V-m^, 

' ° ^ loo 27r y (27r)3 [^2 + (p2 + ;,)2p2 + ^6] 2 ' y > 

The integration over uj leads to a logarithmically-divergent integral over p: 



g2 M p^^p^^Xfdp 



[p2(p2 + ^)2+^6]3/2 

where A is the cut off in the 3-dimensional p space. Although eq. fl2T]) appar- 
ently contains an IR divergence for m/ = 0, this divergence is actually not 
present in this case, if A 7^ 0, since we have then 

< = m^ + ^lnfl + 4^V (22) 



27r2 V A 
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and A plays the role of IR cut off. 

In what follows, the bare mass mg will be chosen such that the renormal- 
ized mass rrib is finite and fixed. This renormalized mass will play the role of 
IR cut off for the calculation of the fermion self energy. 



3.3 Fermion sector and self-consistent equations 



The fermion self energy is calculated from the bare propagator f[T5|) and the 
dressed propagator which is assumed to have the form ([T6|l . such that 

Sy(k) = -A(k-7) -mj. (23) 

Furthermore, if we assume that the dressed scalar propagator has the form 
®,Qb = Gbi where mj, is the renormalized, finite scalar mass fl2T]) . the right- 
hand side of the Schwinger- Dyson equation f|T8|) is (for vanishing frequency 
and after a Wick rotation) 

^ ^ 27r y (27r)3 a;2 + (p2 + ^)2p2 + ^6 



1 

^2 + (p ~k)6 + ml 



X ... , u^B , ^6 - (24) 



This is a convergent integral, and, together with the self energy (l23l) . leads 
to the self consistent equations which must be satisfied by A and m^: 
(i) The equation the fermion dynamical mass should satisfy is obtained by 
taking the trace of the Schwinger-Dyson equation flTSj) . for k = 0: 



{2ixY "^"^ J [co^ + p6 + ml][u^ + (p2 + A)2p2 ' ™6i • (25) 



m 



/J 



If 771/ 7^ 0, the integration over u shows that the dynamical mass must satisfy 

47r2 p p^dp 



AhAfiAh + AfY 



(26) 



where 



Ab = \ P^ + m[ 



Af = ^p^{p^ + \y + mj. (27) 
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(ii) The equation for the coefficient A is obtained by expanding the self energy 
in k, and keeping the hnear contribution only in order to identify it with 
the corresponding term in eq. fl23|l . Using the following equality, valid for any 
function /, 

d^p(k .p){p. 7)/(p^) = -^(k • 7) dp p^+V(/), (28) 

where fio is the solid angle in dimension D, and identifying the coefficients of 
(k-7) in the Schwinger Dyson equation, we obtain the following self consistent 
equation for A 

X = 9^ r ri r^v+AM^ 

where Af,Ab are given in eq. fl271) . Finally, we are left with the two self- 
consistent coupled equations f l26|29l) . which have to be solved simultaneously 
to find the parameters (m/. A) which can be generated dynamically. 



3.4 Numerical analysis and discussion 

In this section, we present our numerical analysis, we comment on the phys- 
ical relevance of our model and the limits of our approximation. 

In order to solve the system of eqs.f l25|) and fl29|) we have to distinguish 
two cases: a) = 0, and h) rrif ^ 0. 

a) If = eq.f l2^ is satisfied automatically, then we have checked 
numerically that eq. fl2^ can be solved with respect to A, for all the range of 
the free parameters g and m^. However, this class of solutions with rrif = 
0, A 7^ is not accepted because of the IR divergences which arise when we 
go to second order approximation in momentum, as we explain in Appendix 
C. For this reason this class of solutions is not presented here. 

b) If we assume that ruf ^ we can divide eq. fl25p by the factor 

to obtain eq.f l26|) in the previous section. Then we can rescale the other 
parameters of the theory with the renormalized mass of the scalar field rrib, 
to obtain the following dimensionless parameters: 

/i = — - I = — 5- e = o , (30) 

nib rrib 2nm'^ 



10 



and the set of coupled equations to solve is, from eqs. 



AhAf{Ah + A 



f) 



+ (31) 



I Jo AlAfiAt + Afj 



where 



Ah = yiT^ Af = v/yu6 + x2(a;2 + /)2. (32) 

We solve the above algebraic system of equations numerically, and a unique 
solution for the pair [l, n) is obtained, if the dimensionless coupling e is larger 
than the threshold Sc — 1.3263. The results for the parameters /^^^ and /i as 
a function of the dimensionless coupling e are presented in fig. ([2]). We would 
like to stress that the singular point = 0, y/l = 0.529 (for e = Ec) is not 
a solution of the system of equations (13T|) because of the restriction /i 7^ 0. 
However the system of equations (13T|) has solutions arbitrarily close to the 
singular point if e > Ec, while for £ < we have no real solutions. 

According to the above results Lorentz symmetry arises for the specific 
model we examine, in the IR limit when p « A^^^. Indeed, from the disper- 
sion relation for the free fermion 



and for p « A^^^, we obtain 

Z?2 ^ ^2 I 



~ p2 ^ ml, (34) 



where we define the rescaled parameters E = u/X and rnf = m^/X that 
correspond to the fermion energy and mass with the correct dimensions [E] = 
[ifif] = 1. 

However, if mj ~ A^^^, the limit p « A^^^ implies that p « mj, and 
the behavior of the particle is then nonrelativistic, the kinetic energy of the 
fermion is given by p^/2mj (note that rrif = fhf for A^/^ ~ mj). We observe 
in fig. ([2]) that there is a small region for which we obtain a relativistic fermion, 
in particular for e > Ec when e is close to the critical value Ec, l^^"^ becomes 
significantly larger than yU . For e » Ec, the mass of the fermion increases 
and becomes comparable to A^^^, this means that the relativistic behavior 
for fermions is restricted to a narrow set of values for the coupling. 
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Figure 2: The parameters ^ = mf/mh, l^/"^ = X^/"^ /mij as a function of e = 
g/{2TTml). The system of equations ([3T]) has a unique solution for e > Sc, which 
for tends asymptotically to /i = and /""^/^ ~ 0.529. For e < we have 

checked numerically that the system pip has no solution. 

The non-perturbative analysis with SD equations, in this article, was 
based on the following two approximations: 1) ladder approximation, 2) first 
order approximation in momentum. 

1) In the case of ladder approximation we make the assumption that the 
bare coupling g is almost equal to the dressed coupling 6. It seems that 
this consideration is consistent with our results, hence we do not expect 
significant corrections in fig. ([2]). In particular, in Appendix B we compute 
the one loop vertex diagram and show that the dressed coupling B*^^^ receives 
small corrections and is therefore close to the bare coupling g. 

2) For 6 ^ Ec, our approximation of taking into account only the first 
order in momentum in the propagators might not be reliable: higher order in 
momentum terms become significantly strong due to IR divergences for ruf = 
0, hence different Ansatze for the propagators should be considered, with a 
larger number of unknown parameters. This would generate a considerably 
more involved numerical problem to solve and is beyond the scope of this 
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article. However, as we explain in Appendix C, we believe that the effects of 
higher orders would smoothen the singular behaviour of the solution when 
e — )■ Eci and that the present singularity is rather an artifact arising from the 
first order approximation in momentum. 

Finally we would like to note that the issue of Lorentz symmetry restora- 
tion becomes problematic when we consider more realistic models with dif- 
ferent species of fermions. In particular, the authors of Ref.|23] find that, 
although renormalization group flows reduce the difference in the velocities 
of different species towards the IR, this is not enough to achieve the exper- 
imental precision of the speed of light, unless we assume an unnatural fine- 
tuning between the Lorentz violating terms of the UV action of the model. 
These results are demonstrated in the framework of perturbation theory by 
considering several simple Lifshitz-type models, and we emphasize that the 
non-perturbative mechanism we propose in this paper suffers from similar 
problems, as we discuss further in the main part of this paper. For example 
if we had considered an extension of our model with two species of fermions, 
with different Yukawa couplings gi and g2, we would have encountered the 
generation of two first order kinetic terms with different Ai and A2 in the 
IR. This would mean that the limiting velocities of fermions would be sig- 
nificantly different in the IR, unless we had chosen gi extremely close to 5^2- 
Such fine-tuning is not consistent with the logic of standard model where the 
Yukawa couplings are significantly different for different species of fermions. 

4 Conclusions 

We considered a 3+1 dimensional model with a Lifshitz-type fixed point 
{z = 3) in which a scalar and a fermion field interact via a Yukawa term. 
The effect of dynamical mass generation, as well as the restoration of Lorentz 
symmetry in the IR limit was examined in the framework of Schwinger-Dyson 
equations. 

An interesting point in this model is that the interaction is super renor- 
malizable, and the only UV divergence present comes from the scalar self- 
energy diagram. Note that, in contrast to the standard case {z = 1) in which 
the divergence is quadratic, the divergence in our model is logarithmic, due 
to the higher powers of momentum in the propagators. In order to absorb 
the UV divergence in our model, we introduce a bare mass for the scalar 
such that our effective theory does not depend on the cutoff of the theory. 
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On the other hand, a bare fermion mass is not necessary in the action since 
a dynamical mass is generated quantum mechanically and is finite. 

Note that the absence of quadratic divergences in the scalar self energy 
diagram sets the hierarchy problem on a new basis, as the scalar field mass 
fiows logarithmically with the UV cutoff, see eq. fl2T]) . However, the absence 
of quadratic divergences in our model can not be considered as a resolution 
to the hierarchy problem as we do not consider all the degrees of freedom of 
the Standard Model (see also [26]). 

The ansatz for the scalar and fermion self energies is based on a linear 
approximation in the external momentum k, and we do not discuss here 
the possibility of generating a Lorentz-invariant kinetic term for the scalar 
field as this term would be of order k"^. The equations arising from the 
Schwinger- Dyson approach are solved numerically and the results are pre- 
sented in fig. ([2]). We find that there is a critical value gc for the Yukawa 
coupling, above which Lorentz symmetry is restored and a mass is generated 
in the low energy limit of the fermionic sector. 

Note, that there are other marginal and relevant operators which can 
be included in the UV action, for example the interaction terms gicj)^^"^, 
^2(f'\l/)^, gscp^, etc . However, had we included these terms in the bare 
action, the corresponding Schwinger- Dyson system of equations would have 
become exceedingly complicated and would in any case have lead us beyond 
the scope of this article. As a first step in a non-perturbative approach, 
we restricted our analysis to the most economical model, namely a Yukawa 
interaction, which is enough to demonstrate the dynamical generation of a 
mass and a first order kinetic term for the fermion. 

Finally, this non-perturbative mechanism for the restoration of Lorentz 
symmetry in models defined at a Lifshitz point may be useful for the study 
of other theories with immediate phenomenological interest, such as QED or 
Higgs models, which are proposed for future investigation. 
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Appendix A: Schwinger-Dyson equation 

The partition function of the theory corresponding to the bare action (ITTI) is 

^[jyVyV] = / ^^[0, V'] 6Xp < iS* + 2 / dtd^X + f]'ljj + ijjT]) 



= exp{iW[j,r],r]]}, (35) 

where j, fj, r] are the sources for (p,ip,ip respectively, and W is the connected 
graphs generator functional. The functional derivatives of the latter define 
the classical fields (pcipc i^c 

— <(h>=(hr. 



6j Z 
6W 1 



6r] Z 
SW 1 



< 7/- >= 



Sr, Z < * 

where 

<••■>= / V[^,ij, ■■)exp\iS + i / dtd^x {jcp + tJV' + [ • (37) 



The proper graphs generator functional r[0c, "^c, "^c] is defined as the Legen- 
dre transform of W, 

T = W - J dtd^x {jcpc + #c + i^cV) , (38) 

where the sources have to be understood as functionals of the classical fields. 
It is easy to check that 

sr 



6T 
6T 



Sip 

s'^T _ /my' 



-J 

V 

-V 

(39) 



15 



The first step for tlie derivation of a self consistent equation involving the 
dressed propagators and vertex is to note that the functional integral of a 
functional derivative vanishes, such that 

|^,)^0. (40) 

Using the different derivatives (l39i) . we obtain then 

^ = + (-A)^ {^^'^,)) - I < 0^ > . (41) 

The vertex, the bare and dressed fermion propagators are respectively 

e ■ ''^ 



where the index refers to vanishing fields, such that a functional derivative 
of eq. fl41l) gives for the fermion self energy 

E,=a7'-G- = -|(A<^^.>) . (43) 







We then express < (pip > in terms of derivatives of W: 



iMc + - < >, (44) 



5j57] Z 



such that 

5 



<0^>) (45) 
6 6^W 



. / 6^W 6r]\ 
\6ri6j5fJ 5j Jq 
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[jt] = -'Gf^Sb, (46) 





and the Schwinger-Dyson equation for the fermion self energy is finally, from 
eq.gSD, 

j:f = tggfegt. (47) 

We stress that this equation is exact, and represents a resummation of all 
quantum corrections, since it involves the dressed quantities on both sides of 
the equation. 



Appendix B: Validity of the ladder approxima- 
tion 

We check here that the ladder approximation is consistent and calculate 
corrections to the coupling constant. This calculation is one-loop like, but it 
takes into account the fermion mass generated dynamically: this is a similar 
approach to Schwinger Dyson equations, and avoids IR divergences. 
This correction is given by the following three-point graph, for vanishing 
incoming momentum: 

3 f du d^p /a;7°-p2(p-7) 



e(i) =g + ighi j 



2tI (27r)3 I ^2 _ p6 _ ^6 _^ ( ^6 _ p6 _ ^6 _^ ' 



and, after a Wick rotation. 



QW =g+^ [ diOp^dp 
TT-^ J 



In the phase where the fermion dynamical mass is generated, mj c::^ mi, (near 
the singular point of fig. |2]but not very close to it), such that 



G^^) ~ g + ^ f duj p^dp 

TT-^ J 



g^ f , 2 J —p^ + ml 
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b 



AtT^ J P "^Py (p6 + ^6)5/2 + ^6)3/2 



The last integral vanishes, and O^^^ ~ As a consequence, corrections to 
the coupling constant can be neglected, and the ladder approximation that 
we used is justified. 



Appendix C: Beyond first order approximation 
in momentum 

In order to go beyond the first order approximation in momentum, we have to 
expand the self energy diagrams up to second order in external momentum. 

In this way we introduce two more parameters Zf and A;,, and we have to 
solve numerically a system of four coupled equations with four unknown 
parameters: 

l_i^mf/mb lf = Xf/ml l^^\^/ml zj = Zf/rub, 

which would satisfy equations of the form 

Fi(ii,lf,lb,Zf,g/ml) = II 

F2{n,lf,lb,Zf,g/ml) = If 

F3{ii,lf,lb,Zf,g/ml) = lb 

Fi{ii,lf,lb,Zf,g/ml) = Zf 

where the functions ^1,^2,^3,^4 are integrals over the energy u and the 
momentum p. Note that these functions could be obtained by expanding the 
fermion and scalar self energies up to second order in momentum k, but the 
corresponding set of equations would lead to a much more involved numer- 
ical problem, and would not change the essential point, i.e. the dynamical 
generation of Lorentz-invariant terms. 
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Nevertheless, one can predict the effect these corrections would have, in 
the vicinity of the critical point of fig|2J In our first order approximation 
in momentum, we disregarded the solutions (/i = 0,// 7^ 0), whereas they 
satisfy our set of self-consistent equations, since // acts as an IR regulator 
when rrif = 0. The reason we disregard these solutions is the following: 
we checked that at nif — )■ 0, the derivatives cPT,f(k)/dk'^ and cPT,h{k) / dk'^ 
become infinite whereas they should be finite, since they are related to wave 
function renormalizations. As a consequence, we believe that the existence of 
this singularity is an artifact of the first order approximation in momentum 
and we expect that second order corrections would smoothen the singular 
behaviour of the dynamical mass. 
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